We have analyzed in detail the renormalization of the axial charge in nuclei
by evaluating the matrix elements of the axial charge operator with
relativistic wave functions, solutions of the Dirac equation with the
relativistic potential, and with non relativistic wave function, solutions of
the Schrödinger equation with an equivalent non relativistic potential. We
have found renormalization effects due to the use of the relativistic wave
functions, enhancing the axial charge in the direction found in earlier
perturbative approaches for nuclear matter.
However, the quantitative results differ from the estimates derived
for nuclear matter significantly. Using the G-matrix derived from a
realistic meson exchange model of the NN interaction [#rupr##1###] a
perturbative estimate of the heavy meson exchange current contribution
to the axial charge at nuclear matter saturation density [#7a##1###] would
yield a renormalization factor of 1.4 and a non perturbative treatment
would lead to enhancement as large as 1.7. For finite nuclei the
enhancement factors considerably smaller, of the order of 1.2 - 1.3.
We argue that this reduction of the renormalization effect is due to
the smaller densities at the surface of finite nuclei, which are
relevant for the evaluation of actual matrix elements. From these
considerations we can also understand the dependence of the
renormalization factor on the momentum transfer and on the transition
actually considered.
The amount of axial charge renormalization depends on the model for the
NN interaction. We have employed a relativistic meson exchange model
(Potential version A of the Bonn potential [#rupr##1###]), which has
been derived to reproduce NN scattering data.
It is fair to quote at this point that using this potential in the
present case there is the assumption that the
relativistic potential constructed to reproduce NN scattering of on shell
nucleons can be extrapolated to deal with negative energy states and on shell
and off shell conditions. This is certainly a strong assumption from which all
the microscopically constructed relativistic potentials suffer, and indeed
different parametrizations of the NN amplitude on shell lead to different
relativistic potentials [#11##1###]. Some efforts have been done to constrain
the relativistic potential to be consistent with the #tex2html_wrap_inline1414#N elementary
amplitudes [#12##1###] and this leads to potentials like the one obtained here
but about one half their strength. Even then this potential is constructed at
the level of the impulse approximation or low density limit, tρ, and
many body effects should modify it. It is clear that many efforts are still
necessary to be able to claim that an unambiguous microscopical relativistic
potential has been determined. On the other hand one can take a more
phenomenological approach and say that a certain relativistic potential has a
wide degree of phenomenological success, providing fair nuclear binding
energies, spin-orbit splitting, nucleon nucleus cross sections and
polarization observables, etc. [#10##1###,#klein##1###].
The potential we have used is one of such and provides
empirical support for the axial charge renormalization found,
but this does not
exclude the possibility of other potentials with the same degree of
phenomenological success and still providing different axial charge
renormalization. The ultimate answer to this question is tied to the progress
in our understanding of the meaning and accurate strength of the relativistic
potential. Meanwhile, by using a fair and plausible model we have done
detailed calculations and showed that the results are sufficiently different
from the perturbative results to encourage the use of the present approach in
future works dealing with the problem.
Two of us, A. Gil and E. Oset wish to acknowledge the hospitality of the
University of Tübingen and H. Müther the one of the University of Valencia.
E. Oset acknowledges support from the Humboldt Foundation. The work has been
partially supported by the EU, program, Human Capital and Mobility contract
no. CHRX-CT 93-0323, the CICYT contract no. AEN 93-1205
and the Graduiertenkolleg ``Struktur und Wechselwirkung von Hadronen
und Kernen'' of the Deutsche Forschungsgemeinschaft (DFG Mu 705/3)
<#370#>Appendix: Matrix elements of the axial charge operator.<#370#>
<#371#>A) Relativistic case:<#371#>
We write here the matrix element for the #math57#γ0γs operator between relativistic wave functions
#math58#
;SPMlt; n'l'j'm'| γ0γ5ei#tex2html_wrap_indisplay1418##tex2html_wrap_indisplay1419#| nljm ;SPMgt; #tex2html_wrap_indisplay1420#(a.1)
We distinguish 4 cases
#math59#
a)j' = l' + 1/2, |
j = l + 1/2 |
b)j' = l' + 1/2, |
j = l - 1/2 |
c)j' = l - 1/2, |
j = l + 1/2 |
d )j' = l - 1/2, |
j = l - 1/2 |
and the resulting matrix element is
#math60#
#tex2html_wrap_indisplay1431##tex2html_wrap_indisplay1432#iλ(- i)#tex2html_wrap_indisplay1433#r2dr#tex2html_wrap_indisplay1434##tex2html_wrap_indisplay1435##tex2html_wrap_indisplay1436# - #tex2html_wrap_indisplay1437##tex2html_wrap_indisplay1438##tex2html_wrap_indisplay1439#jλ(qr)
#math61#
(2λ +1)1/2Yλ, m'-m*(#tex2html_wrap_indisplay1441#)Ai#tex2html_wrap_indisplay1442#(a.2)
where Ai is given for each of the cases a) b) c) d) listed above by
#math62#
Aa = #tex2html_wrap_indisplay1445#C(j + 1/2, λ, j' - 1/2;000)
#math63#
{(j' + m')1/2(j + 1 - m)1/2C(j + 1/2, λ, j' - 1/2;m - 1/2, m' - m)
#math64#
- (j' - m')1/2(j + 1 + m)1/2C(j + 1/2, λ, j' - 1/2;m + 1/2, m' - m)}#tex2html_wrap_indisplay1448#(a.3)
#math65#
Ab = #tex2html_wrap_indisplay1450#C(j - 1/2, λ, j' - 1/2;000)
#math66#
{(j' + m')1/2(j + m)1/2C(j - 1/2, λ, j' - 1/2;m - 1/2, m' - m)
#math67#
+ (j' - m')1/2(j - m)1/2C(j - 1/2, λ, j' - 1/2;m + 1/2, m' - m)#tex2html_wrap_indisplay1453#(a.4)
#math68#
Ac = #tex2html_wrap_indisplay1455#C(j + 1/2, λ, j' + 1/2;000)
#math69#
{(j' + 1 - m')1/2(j + 1 - m)1/2C(j + 1/2, λ, j' + 1/2;m - 1/2, m' - m)
#math70#
+ (j' + 1 + m')1/2(j + 1 + m)1/2C(j + 1/2, λ, j' + 1/2;m + 1/2, m' - m)}#tex2html_wrap_indisplay1458#(a.5)
#math71#
Ad = #tex2html_wrap_indisplay1460#C(j - 1/2, λ, j' + 1/2;000)
#math72#
{(j' + 1 - m')1/2(j + m)1/2C(j - 1/2, λ, j' + 1/2;m - 1/2, m' - m)
#math73#
- (j' + 1 + m')1/2(j - m)1/2C(j - 1/2, λ, j' + 1/2;m + 1/2, m' - m)}#tex2html_wrap_indisplay1463#(a.6)
<#434#>B) Non relativistic case:<#434#>
we evaluate matrix elements of the #math74##tex2html_wrap_inline1465#(#tex2html_wrap_inline1466# + #tex2html_wrap_inline1467#)/2M operator between non relativistic states
#math75#
;SPMlt; n'l'j'm'|#tex2html_wrap_indisplay1469#(#tex2html_wrap_indisplay1470# + #tex2html_wrap_indisplay1471# ')/2M| nljm ;SPMgt; #tex2html_wrap_indisplay1472#(a.7)
The derivation of this matrix elements requires a bit more algebra than the
non relativistic case. With the help of some useful formulas from the
appendix of ref. [#13##1###] we obtain the following result
#math76#
;SPMlt; n'j'l'm'|#tex2html_wrap_indisplay1474##tex2html_wrap_indisplay1475#ei#tex2html_wrap_indisplay1476# #tex2html_wrap_indisplay1477#| nljm ;SPMgt; =
#math77#
i(- 1)j'+l'+1/2#tex2html_wrap_indisplay1479##tex2html_wrap_indisplay1480##tex2html_wrap_indisplay1481#C(jλj';1/2, 0, 1/2)
#math78#
C(jλj';m, m' - m)iλYλ, m'-m*(#tex2html_wrap_indisplay1483#)#tex2html_wrap_indisplay1484##tex2html_wrap_indisplay1485##tex2html_wrap_indisplay1486#
#math79#
F(n'l'j', nlj;λ)#tex2html_wrap_indisplay1488#(a.8)
with #math80#λ + l + l' an odd number, where the last function is given by
#math81#
F(n'l'j', nlj;λ) =
#math82#
δj, l+1/2(l + 1)1/2#tex2html_wrap_indisplay1492#r2drφl'(r)[#tex2html_wrap_indisplay1493# - #tex2html_wrap_indisplay1494#φl(r)]jλ(qr)
#math83#
- (- 1)j-j'+l-l'δj', l'+1/2(l' + 1)1/2#tex2html_wrap_indisplay1496##tex2html_wrap_indisplay1497##tex2html_wrap_indisplay1498#
#math84#
#tex2html_wrap_indisplay1500#r2dr[#tex2html_wrap_indisplay1501# - #tex2html_wrap_indisplay1502#φl'(r)]φl(r)jλ(qr)
#math85#
- δj, l-1/2l1/2#tex2html_wrap_indisplay1504#r2drφl'(r)[#tex2html_wrap_indisplay1505# + #tex2html_wrap_indisplay1506#φl(r)]jλ(qr)
#math86#
+ (- 1)j-j'+l-l'δj', l'-1/2l'1/2#tex2html_wrap_indisplay1508##tex2html_wrap_indisplay1509##tex2html_wrap_indisplay1510#
#math87#
#tex2html_wrap_indisplay1512#r2dr[#tex2html_wrap_indisplay1513# + #tex2html_wrap_indisplay1514#φl'(r)]φl(r)jλ(qr)#tex2html_wrap_indisplay1515#(a.9)